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An effective medium approach together with a multiple scattering formalism is considered to 
study the steady-state dynamics of suspensions of spheres and rigid stiff polymer chains without 
excluded volume interactions. The polymer chains are taken to be so long that Gaussian statistics is 
applicable. We have considered the dynamics of probe objects in a solution containing spheres and 
polymers. The probe object is either a sphere or a polymer. We have studied different conditions 
of the solution where some or all of the spheres are frozen in space. The effective medium equations 
have been solved self-consistently for finite volume fractions of spheres $ SP and polymers <E> POi , 
respectively, and the important dimensionless variables that are controlling the dynamical behavior 
have been identified. In particular,the role of the geometrical parameter t = — (o is the radius 
of any sphere and R g the radius of gyration of a polymer chain) is discussed. The translational 
diffusion coefficients of the moving probe sphere D s and of the probe polymer chain D p ', and 
the shear viscosity of the suspensions have been derived. When the polymers are present in the 
solution, both the friction coefficient of the labeled mobile sphere or that of the probe polymer chain 



and the shear viscosity diverge as $ 



POL 



0.31. Also, when polymers are diffusing in a suspension 



of fixed spheres an optimum range of <E> SP that maximizes the difference in the diffusion coefficients 
of polymer chains characterized by distinct t values has been noticed. 



I. INTRODUCTION 



Starting with the pioneering work of Einstein |jj , a lot 
of effort was devoted to developing the hydrodynamics of 
suspensions at finite concentrations, owing to both the in- 
trinsic challenge of the problem and to the considerable 
utility of such a theory in the rheology of concentrated 
polymer suspensions. Most of the work done is related to 
the hydrodynamics of suspensions of hard spheres, specif- 
ically to the computation of virial-type expansions — in 
the volume fraction <!> (low but non-zero) of the sus- 
pended particles — for the change in viscosity Sq and in the 
friction coefficients ( (translational and rotational). An 
incomplete review of the development of this field could 
include the divergence-free theories of Peterson and Fix- 
man H and Batchelor and Green |i| , the multiple scatter- 
ing method elaborated by Freed and Muthukumar ||||] 
and Muthukumar and Freed the expansion in cor- 

relation functions used by Beenaker Q and the angular 
momentum diagrammatic expansion combined with the 
multiple scattering technique in Thomas and Muthuku- 
mar M. 

The viscosity and frictional properties of infinitely di- 
luted polymer solutions received a wide interest (see Ref. 
ju|), but few attempts were made to construct a detailed 
theory for finite concentrations |T^-jl^] , and even less the- 
oretical results are available for dispersions of distinct 



types of elements. 

In this paper we derive explicit expressions for the 
change in viscosity due to suspended objects and for 
the translational diffusion coefficients, as functions of the 
solid filling fractions of spheres and of polymer chains, 
<J> 5P and $ POL , respectively, for a suspension of poly- 
mers and rigid, fixed spheres, in the stationary, time- 
independent limit. Qualitative features of the hydrody- 
namic behavior of the mixed suspensions are discussed 
and it is shown that when the polymers are present in 
the solution, the dynamics is frozen at $ POL ~ 0.31. 
Results in the limits of very low and very high poly- 
mer concentrations (Rouse) and also for the dynamics of 
a random array of spheres m| , are recovered. We have 
found that the interaction between polymers, spheres and 
the fluid is controlled by the dimensionless parameters 
/3= x = | and t = Zl, where f = £(<P SP , $ POi , t) 
is the hydrodynamic screening length. 

Technically, our calculations are based on the multi- 
ple scattering formalism P,p|,p^|J^] and on the assumption 
that the physical properties of the suspension, averaged 
over the random position distribution of the particles, 
are equivalent to those of an homogeneous effective fluid 
characterized by the effective viscosity r] eS and by the al- 
tered Oseen propagator of the velocity field Q(r, r') (to be 
defined later). 

The main simplifying hypothesis is the absence of 
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excluded volume interactions — sphere-sphere, polymer- 
polymer, or polymer-sphere — both the polymer chains 
and the spheres being penetrable, ghost-type objects. 
Also, the solvent is assumed incompressible and de- 
scribed by a linearized Navier-Stokes (N-S) equation. 

Other assumptions made and the basic features of the 
theory are : 

a) the suspension of spheres and polymers is monodis- 
perse, all spheres having the radius a and all polymers 
having the length L. 

b) except the sphere producing the flow, all others are 
fixed. All polymers are Gaussian. 

c) the specific calculations are valid only in the 
limit of steady-state shear flow in the linear response 
regime, when the average force (F(r)) exerted upon the 
fluid is given by : 

(F(r)) = (K*(v))(r) 

where K is a tensor-like quantity and * is the convolution 
operator: f * g = J f(r — r')g(r') dr'. 

d) both the spheres and the polymer chains are ran- 
domly distributed, and after the configurational average 
is taken the suspension becomes isotropic. 

e) we employ the usual preaveraging approximation 
jlO| , where one replaces the value of the effective Oseen 
tensor Q(r,r') and of its generalized inverse with their 
configurational averages (<7(r, r')) and (i/ _1 (r, r')\. 

f ) two models have been analyzed : 

S-P an ensemble of rigid spheres in a dispersion of polymers. 
All but one sphere are fixed. The tracer sphere gener- 
ates the flow. 

P-S a collection of polymer chains immersed in a suspension 
of immobile spheres, with an external flow imposed. 

The hydrodynamic interaction between fluid and par- 
ticles is implemented by enforcing stick boundary con- 
ditions. In each example we calculate the translational 
diffusion coefficient of the test object ( a polymer chain or 
the tracer sphere ) and the change in the shear viscosity 
of the solution due to the presence of the added particles, 
as functions of & SP , $ POi and the coupling parameter 
t = ^2.. 

a 

The paper is organized as follows : next section deals 
with the formal theory of the steady-state velocity flow in 
a solution of Np polymer chains and fixed spheres , in 
Section III we apply the derivation to the separate cases 
of polymers and spheres, which are needed in Sec. |y| 
where we solve and show the results for the four specific 
examples mentioned previously. Finally, in Section 
we present, discuss, and review the main findings of our 
paper. For the reader interested in reproducing and ex- 
tending our results, detailed calculations are displayed in 
the Appendices. 



II. FORMAL THEORY OF THE VELOCITY 
FIELD IN A SUSPENSION OF POLYMERS AND 
SPHERES 



To describe the dynamics of a hydrodynamic system, 
a knowledge of the velocity field v(r,t) produced by the 
forces acting upon the fluid is required. Our goal is to find 
an expression for v(r) (stationary limit) for the system 
we are interested in and to extract the appropriate phys- 
ical information (diffusion coefficient, viscosity). This is 
accomplished using the multiple scattering technique as 
shown below. 



A. Fluid and Particles 



Under the assumptions indicated in the previous sec- 
tion and taking the fluid density equal to unity, the sta- 
tionary velocity field for an incompressible, newtonian 
viscous fluid is governed by the linearized N-S equation : 



?7oAv(r) + Vp(r) = F(r) 



(1) 



where F(r,t) and p(r, t) are the external force driving the 
flow ( including the random thermal contribution ) and 
the pressure, and 770 is the kinematic shear viscosity coef- 
ficient. Introducing the Fourier transform of any function 
A(r) as : 

A(r) = J 7^~vj A(k) exp(— ik ■ r) = J A(k) exp (— ik ■ r) 
A(k) = /dr A(r)exp(ik- r) (2) 

one finds the formal solution of Eq. (|l|) as a convolution : 
v(r) = (G * F)(r) (3) 

with G(r,r') being the Oseen tensor jl6| (in dyadic no- 
tation) : 



GOO = ^ 



87r?7o|r — r'| 



Here, 1 is the unit tensor and k is the unit vector pointing 
in the k direction. 

Physically, G(r, r') — the force propagator — is convey- 
ing the hydrodynamic disturbances, converting a point- 
like force acting at r' to the velocity of a fluid particle 
located at r. When a fixed object is introduced in the 
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fluid, the velocity flow is perturbed due to the scatter- 
ing off the particle surface. One needs then to include in 
Eq. ([!]) a term ~p OBJ ^ FL (r) corresponding to the force 
exerted by the particle upon the fluid : 



r/ Av(r) + Vp(r) = F(r) 



? OBJ^FL 



(r) (5) 



To couple in a simple way the dynamics of the solvent and 
of the object, we choose no-slip boundary conditions : 



FLUID 



v (r s ), r s 6 object surface. (6) 



In Eq. (0) F(r) is zero inside and on the surface of 
the particles and ~f obj ~* fl {y) is non-zero only on the 
surface of the suspended particles. The equation of mo- 
tion is valid throughout the volume of the fluid. Because 
we want to replace the fluid and the suspended particles 
with an effective medium characterized by an average ve- 
locity field defined everywhere, it is necessary to apply 
Eq. ([i]) inside the volume of the particles as well. This 
is true provided the full stress tensor II(r) |lj],|| obeys ( 
neglecting the inertial term in the equation of motion of 
the object ) : 



V • II(r) = 
V • II(r s ) = F OBJ ^ FL (r s ) 



, r 6 inside object. 
, r s 6 object surface. 



(7) 



The formal solution of Eq. (||) can be written as : 

v(r) = G*F| r + G*F OBJ ^ FL \ r (8) 

Using the boundary condition (Q), we can eliminate the 
unknown quantity ~f? OBJ ^ FL (r) — as shown explicitly in 
the next section — and arrive at the exact equation that 
gives the microscopic velocity field everywhere in the 
fluid : 

v(r) = G*F| r + G*T OBJ *G*F| r (9) 
where the multiple convolution product means : 
f * g * h| r = Jdr dr' f (r - r') • g(r' - r") • h(r") 

and T OBJ represents the flow propagator, a tensor op- 
erator that transforms the velocity field incident on the 
object in a force field located on its surface and acting 
upon the fluid. T OBJ depends in general on the position, 
structure and geometry of the object. The interpretation 
of (^) is straightforward : there are two contributions to 
the velocity of the fluid v(r) at point r, a direct wave due 
to the external force F acting at some point r', and an 
indirect contribution coming from the scattering off the 
surface of the object of the disturbance produced by F 
at r' (see Fig. 1). 

When many objects are present in the solution, one 
should add all the contributions coming from all possible 



scattering sequences but, even for two bodies only, Eq. 
(0) becomes a series with an infinite number of terms, 
the convergence of which might be problematic. For 
some specific cases (polymers, spheres) these dif- 

ficulties can be overcome, which was the reason we de- 
cided to investigate solutions of mixed polymer chains 
and spheres. For example, when two objects labeled a 
and b are present in the solution, the multiple scattering 
expression for the velocity field v(r) will be : 

v(r) = G * F| r + G * {T a + T b } * G * F| r 

+ G*{T a *G*T b + T 6 *G*T a }*G*F| r (10) 

+ G * {T a * G * Tt, * G * T Q 

+ T b *G*T a *G* T b } * G * F| + . . . 



Being interested only in the average properties of the 
system, we will average in Eq. (|^) upon the random posi- 
tion of the body, thus obtaining the average velocity field 
u(r) : 

u(r) = (v(r)) =i |dR°v(r)= G * F| r 

+ G*(T aBJ )*G*F| r (11) 

with R° denoting the position of the center of mass of 
the object. Note that if excluded volume interactions are 
considered, appropriate particle distribution functions 
should be used. 

To make the connection with the experimentally mea- 
surable quantities, we will turn to Eq. (^J) and average it 
directly over the particle distribution to get the Navier- 
Stokes equation describing the effective fluid : 



7 ?0 Au(r) + V (p(r)> - £ * u| r = F(r) 



(12) 



In writing the equation above we introduced a new quan- 
tity, the £(r, r') operator, called the self-energy of the 
fluid and defined by the relation : 



£ * u| 



OBJ^FL 



(13) 



Note that this relation implies a linear response regime 
of the shear flow. 

£(r,r') is the essential quantity that encompasses all 
the information regarding the change in the viscoelastic 
properties of the fluid due to the presence of the object(s). 
To illustrate this, let us introduce the Fourier transform 
of the self-energy, knowing that the averaged solution is 
translationally invariant : 



£(k) 



J dr£(r 



r') exp(ik • (r — r')) 



(14) 



Then Eq. (ft) becomes : 
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[rj k 2 t - S(k)] • u(k) - ikp(k) = f(k) 



(15) 



Using the incompressibility condition k • u(k) = and 
decomposing £(r, r') into its transverse and longitudinal 
parts applying the relations : 

S(k) = S_L(k) (1-kk) +Sy(k)kk 

(i - kk) ■ (i - kk) = (i - kk) (i6) 

kk • (1 - kk) = 

we can eliminate the pressure p(k) to solve Eq. ([l5]) for 
the Fourier component u(k) of the averaged velocity field 
u(r) : 



u(k) = 0(k) • F(k) 



%A; 2 -Sx(k) 

We recover then, immediately : 
u(r) = g*F\ r 



(17) 



(18) 



The interpretation of £ emerges clearly. What we 
achieved is to replace the initial system fluid + object(s) 
with an effective medium of a certain viscosity T) where 
the hydrodynamic disturbances are propagated by the 
modified Oseen tensor G(k), with all the complexities of 
the multiple scattering processes involving the suspended 
object(s) captured by the S tensor. 

In the hydrodynamic limit |k| — > 0, the total change 
in viscosity is given by : 



rf 



off 



— = - lim — -J^rY,j_(k) 

rj k^o 770 ok z 



(19) 



in the solution. Note that in general the friction coeffi- 
cient is a tensor-like quantity, in which case the previous 
expression should be adjusted correspondingly. The di- 
mensionless variable that we have computed and plotted 
in this paper is the relative translational diffusion coeffi- 
cient : 



(21) 



with D® the translational diffusion coefficient of a parti- 



cle in the pure fluid (e.g. U( 



k B T 



for a sphere and 



D 



_ 8V2 



k B T 



for a single polymer chain ). In 



* 3 6irV2iiTioRg 

the remaining of this work, we will denote this relative 
diffusion coefficient simply by D. 

Having established that S(k) contains the relevant dy- 
namical information , we need to actually compute it by 
relating it to the known quantity T OBJ . Here lies the 
essence of the multiple scattering formalism. Averaging 
Eq. (||) directly and using the definition of X from Eq. 
( |l3| ) we can rewrite the solution for the averaged velocity 
field u(r) in terms of the G propagator of the pure fluid : 



u(r) = G * F| r + G * E * u|. 

Iterating indefinitely, one obtains : 

u(r) = G * F| r + G * £ * G * F| r 
+ G*£*G*£*G*FI 



(22) 



(23) 



Representing X as a series in the number i of distinct 
scattering events : 



£(r,r') = ^E«(r,r') 



(24) 



If £_i_(k = 0) 7^ , the hydrodynamic interaction is 
screened, with the screening length £ given by : 



r 



-S±(k = 0) 



(20) 



Similarly, formulas for the friction coefficients C (transla- 
tional, rotational, cross translational-rotational ) of the 
object can be derived from S(k). We will explicitly show 
this, for the translational friction coefficient, in the fol- 
lowing sections. 

More practical, from an experimental point of view, 
is the diffusion coefficient D of a particle moving in a 
solution, calculated from the Einstein formula : 



k B T 

cm 



where T = temperature, ks is the Boltzmann constant 
and $ is the volume fraction of the suspended particles 



substituting in Eq. (|23J) and comparing to the formal so- 
lution from Eq. (|ll]), it follows that : 



£«(r,r') = (T OBJ (r,r')) 
£(2)( r , r ') = _ (T OBJ ) * G(r l7 r 2 ) * <T OSJ )| 

S( 3 > (r, r') = (T OBJ ) * G * (T OBJ ) * G * <T OBJ ) 



(r,r') 



(25) 



(r,r') 



where integration over ri and r 2 in the multiple convo- 
lution is understood. 

The expressions displayed above are particular for the 
case of only one particle in the solution, but they illus- 
trate the basics of the formalism used in this paper. Now 
it is straightforward to generalize the procedure devel- 
oped before to deal with a system where more than one 
object, or different types of particles are immersed in a 
fluid. 
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Finally, at any volume fraction $ of the solute par- 
ticles, we compute self-consistently the contribution 
S(r, r') of the suspended particles to the divergence of 
the full stress tensor of the fluid. This is accomplished 
by assuming linear superposition of the contributions of 
each particle to the self-energy of the fluid (assumption 
which is exact when there are no thermodynamic corre- 
lations among particles, as in the present work). Then, 
for a suspension of N particles, one can replace N — 1 
of them with an effective medium described by the mod- 
ified Oseen tensor Q(r^r'), which is a function of some 
Ejv-i (see Eqs. (|lj]r|[ij)). Then, it follows that <5£, 
the self-energy increase due to the remaining particle is : 

8X(g(X N _ 1 )) = ^-± (26) 

The left-hand side of the equation can be evaluated as 
outlined for the case of one object in a fluid, but using 
the Q propagator instead of G of the pure fluid. Thus 
we have obtained a self-consistent equation for calculat- 
ing the full Sjv-i- 



B. Calculation of S(r, r') for a suspension of N5 
spheres and Np polymer chains 

We consider Ng fixed, rigid spheres, each of radius a, 
and Np mobile, Gaussian polymer chains, each of length 
L = nl (I = Kuhn length; n = number of beads ), im- 
mersed in an incompressible fluid of viscosity 770, in the 
regime of stationary flow. 

The static velocity field in the suspension is described 
by the N-S equation together with stick boundary condi- 
tions : 

N P ,n 

-7| Av(r) + Vp(r) = F(r) + ^ <5(r-R m )<x m 

CK,Z>1 

+ dn b 5(r-R b )<T b (n b ) (27) 

b=i ■> 

V • v(r) = 

R ai = u a + Lj a x S OT = v(R ai ) (28) 
R fc = = v(Rb) 

Here, the greek symbols a, f3, . . . are labeling the polymer 
chains and are running from 1 to Np, i, j, k, . . . are indices 
(g {1, n}) for the beads on any arbitrary chain and 6, c . . . 
are labels for the spheres (range 1, N5). R Q ; — R" + S ai 
is the position vector of the i-th bead of the chain a, 
with R^ denoting the position vector of the center of 
mass of the a chain and with S a i the position vector of 
the z-th bead with respect to the center of mass of a. 



R5 = R" + r b (£l b ) is the decomposition of the position 
vector of a point on the surface of the sphere b in the 
position vector of its center of mass and the relative co- 
ordinate of the surface point in the center of mass frame 
(\v b \ = a, for any sphere); Q indicates the orientation of 
the vector v b . cr a i and <r b (£l b ) are the densities of force 
exerted by the i-th bead of chain a and by the sphere 
b at the surface point rf,, respectively, upon the fluid. 
When there is no risk of confusion, we will write only cr b , 
the dependence of Q b being understood. Finally,F(r) is 
some force acting at r that generates the flow (e.g. the ef- 
fect of any sphere or polymer chain moving with uniform 
velocity) . 

From the boundary conditions for the polymer chains 
and because we are neglecting the inertial terms in the 
equations of motion for the polymers, the total force and 
torque acting upon any chain must vanish. We can write 
then : 

n 

- ^ o- ai = 

i=l 

n 

- S «i x Vcci = , a = ITNp" (29) 
i=i 

The formal solution of Eq. (^) is : 

Np,n 

v(r) = G*F| r + ^2 G(r— R m ) ■ a ai 

a,i=l 

N S , 

+ dQ b G(r - R b ) ■ tr b (Sl b ) (30) 
6=1 

Using the boundary conditions ( p8| ) and the constraints 
( p9| ) we can eliminate the unknown forces cr a i and cr b , 
to express the solution for the velocity field as a multiple 
scattering series in terms of the single-object flow prop- 
agator T Q for the chain a and T;, for the b sphere (see 
Appendix A for details) : 

v(r) = G(r) * F| r + ^ G * {T« + T fc } * G * F| r + 

a, b 

+ Y J G *{ T <**G*T b + T b *G*T a }*G*F\ r 

a, 6 

+ G*T a *G*T /3 *G*F| r (31) 
+ G*T fc *G*T c *G* F| r + . . . 

b^c 

where the sum should be continued over all possible scat- 
tering sequences. In this expression, the first factor on 
the right-hand side represents the direct wave, the sec- 
ond contains the single scattering processes from only one 
polymer chain or one sphere, the third includes the se- 
quences sphere-polymer and polymer-sphere, etc... . Any 
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sequence of T operators is valid except those involving 
two consecutive scatterings off the same sphere or the 
same polymer (the exclusion constraint). 

The macroscopic equation for the effective velocity 
field u(r) is retrieved by performing a configurational 
(position) average over the distribution of the particles, 
mathematically expressed as : 



(32) 



with the "0" superscript referring the center of mass and 
< >ai,aj,... being an average over the distribution of the 
segments of the a chain about its center of mass. The 
probability distribution function is taken to be Gaussian 
(see [0). Then u(r) is calculated from : 



u(r)= (v(r)) = G * F + £ G * {(T Q ) + (T 6 )} * G * F 

cx,b 

+ G * {( T " * G * T &) + ( T fe * G * T a )} * G * F 

a, 6 

+ Y G * (T Q * G * T ) * G * F (33) 



+ ^2 G * (T& * G * T c ) * G * F + . . . 

b^c 

where all convolutions are actually functions of r. 

As outlined in the previous section, the meaningful 
physical quantity is the self-energy tensor S(r, r') defined 
as : 

£*u = /^<5(r- R Qi )er OT + ^ J dtt b 5(r -R b )cr b \ 



£ POL *u + S sp *u 



(34) 



Next, the solution of the configurationally averaged 
Navier-Stokes equation is still Eq. (p3|) , but now the to- 
tal X includes the effects of both the polymers and the 
spheres. Expanding once again the self-energy in the 
number of distinct scattering events (i.e. the number of 
T Q .fc operators) and comparing the terms with the same 
number of T factors in Eqs. (^3|,[33]), it follows that : 

£«(r,r') = ]T<T Q )+£<T 6 > 

a b 

£( 2 > (r, r') - ]T < T « * G * T /3> + E < T " * G * T c ) 

a^(3 b^c 

+ ^{(T Q * G * T fc ) + (T(, * G * T Q )} 

- ( T «> * G * ( T /3> - £ ^ * G * < Tc > ( 35 ) 

a,/3 b.c 

- ^{(T Q ) * G * (T b ) + (T b ) * G * (T Q )} 

a . b 

(36) 



Assuming that the spheres do not rotate , we can de- 
rive the following expressions for the T a . b operators (see 
Appendix A) : 



T Q (r,r') =-^<5(r-R Qi ) [K" 1 ^, S aj ) 

n 

- J2 K- 1 (S m ,S Q/ )- g r 1 -K- 1 (S aP ,S QJ )],5(r-R a 



i,i'=i 

+ rotational terms 



(37) 



i,j 

T b (r - r') = - Jdn b dn' b 5{r - R b ) K^(il b , fi' b ) 6(r' - R' 6 ) 



in which R;, and TL' b are the position vectors of separate 
points on the surface of the same sphere labeled b and 
K" 1 and K^" 1 are the generalized inverse operators for 
the single polymer chain and the single sphere |]l2| , |l5| , 
defined by : 



j^K-^Si.Sj-) • -S fc ) = lS ik (38) 

Jdn'i K _1 (f2b, n'l) ■ G(r b (n' b ') - r b (n' b )) = a s(n b - n' b ) 



To avoid an increasingly intricate notation, we will adopt 
the following short-hand notation : 



K (S a ^, S Q j) > K 

K-\n b ,n' b )^K bb , 

G(R a j — Rh) — > G a j :b 



i 

1 



(39) 



The Einstein summation convention over repeated in- 
dices is implied everywhere but where explicitly not fol- 
lowed . We also drop the index a when referring to quan- 
tities not specifically depending on a particular polymer 
chain. 

Some remarks are necessary before embarking upon 
some concrete calculations. In the limit of our static de- 
scription of the suspension of polymer chains and spheres, 
Eqs. ( |35| ) are exact, for noninteracting as well as for in- 
teracting objects, and will describe the dynamics of the 
suspension with any degree of accuracy, although prac- 
tically it may require a strenuous effort. To obtain an- 
alytical results, we are limiting ourselves to the case of 
noninteracting polymers and spheres and, further, we will 
approximate the total self-energy tensor S(r, r') with the 
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first term in the expansion ( |3q ) that contains the contri- 
bution from single scattering events : 



S(r,r')^SW(r,r') 



(40) 



but, for the convenience of notation we will still denote 
it by X(r,r'). Actually, for the noninteracting situation, 
this is the leading term as can be seen from Eq. (|35|) by 
breaking the averages, with the contribution from higher 
order processes vanishing due to the lack of correlations 
among particles. 

As already pointed out, S(k) contains the sought infor- 
mation about the transport properties of the suspension. 
To evaluate the self-energy, first we take the configura- 
tional average (32) over the T operators in Eq. (p7|), then 
we Fourier transform their expressions, obtaining : 

n 

<T Q (k)> = -- £ (exp (ik • (S, - S,)) [K-i QJ 



,V>1 



(41) 



(T b (k)) = - i y"rfr! h dfi fc K^ 1 , exp [ik • (r b - r' b )] (42) 



The preaverage over K 1 part produces terms like 
(Ky )., (exp (ik ■ (Si — Sj)))^. — with the ij indices sig- 
nifying an average over the distribution of the segments 
of any chain — so S(k) can be written as : 

n 

S(k) ~ E«(k) = -c POL J2 (exp (* • (S, - Sj))) y 



(43) 



J,H'>1 

c SP J dQ dtf K7 l (Q,, Q!) exp [ik ■ (r(fi) - r'(fi'))] 



where r(f2) , r'(fi') are two generic points on the surface 
of any sphere and c POL = ^p- , c SP = ^ are the concen- 
trations of the polymers and of the spheres, respectively. 



III. CALCULATION OF S(k) AND EFFECTIVE 
MEDIUM THEORY FOR SOLUTIONS OF 
SPHERES OR POLYMERS ONLY 



In this part of the paper, the problems of constructing 
the effective media in the separate cases of a suspension 
of spheres and a polymer dispersion are addressed. 



A. Calculation of W(k) for Ns spheres immersed in 
a fluid : a review 



This problem is well documented in the literature |[7|Jl 5|] 
so we will quote the main results we need. The system 
consists in N<j, rigid, penetrable spheres of radius a, im- 
mersed in an incompressible, newtonian fluid of viscosity 
t)q. All spheres but one, which is moving at constant ve- 
locity and generates the flow, are fixed. The suspension 
is assumed stationary. 

A remark should be made regarding the notation. Be- 
cause we want to distinguish this particular system from 
the others ,we will denote the self-energy of the fluid in 
the presence of the spheres only by W SP instead of X SP . 

Along the lines previously exposed,we can map the sus- 
pension averaged over the uniform distribution of the 
spheres to an effective fluid of viscosity rf<? p where the 
disturbances are propagated by the modified Oseen ten- 



sor g SP (k) 



Q SP (k) = 



t-kk 



rjok 2 - W[ p (k) 



(44) 



with the transverse part of the self-energy computed from 
Eq. ©. 

In order to solve self-consistently for W SP (k) we ex- 
pand its transverse part, in the small k limit : 



w[ p (k) = (-^s 2 P -mwr(^ SP )k 2 ) 

Then the modified Oseen propagator becomes : 

1-kk 



(45) 



Q SP (k) 



m(i + w? p )(k 2 + ?s P ) 

ssp- 



1 + W( 



SP 



(46) 



where Wf p is some function of & SP . All other symbols 
were already defined. 

W SP (k) has the following form (see Appendix B) : 



W SP (k) = -Girrioac 



SP 



1 



— (l + Wf p )l + 



Ii/ 2 (f)Ki/ 2 (f) 
1 

8 2 

9I 3/2 (f )K 3/2 (f ) " 5l 1/2 (f)K 5/2 (f)_ 



+ ^ SP (l + wr)k 2 



(47) 



%* SP (1 + Wf)k 



SP\ 1,2 



27I 3/2 (f)K 3/2 (f) 



5I 1/2 (fOK 5/2 (#) 
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where I and K are the modified Bessel functions of the 



first kind, <E> 



SP 



z SP is the volume fraction of 



spheres and z is the versor of the Oz direction. Pro- 
jecting out the transverse part of W SP (k) by dotting in 
with (1 — zz) and comparing the result with the expan- 
sion (J45|), one arrives at the following set of equations for 
x = §j and Wf p : 



x 2 = »$SP 



1 



l 1/2 (x)K 1/2 (x) 

w 1 sp (<p sp ) = |i + ^ SP 



1 



_li/ 2 {x)K 1/2 (x) 
2 



9I 3 /2(^)K 3 / 2 (a;) 5I 1/2 (a;)K5 /2 (x) 



(48) 



As stated before, the readily measurable quantity we 
calculated is the relative diffusion coefficient of the mo- 
bile sphere (relative to the bare diffusion constant), in- 
troduced in (|2l[). From (|is| ) we get immediately : 



B SP = D sp t = I 1/2 ( X ($ sp ))K 1/2 (x(4> sp )) 

1 

W ? P (<PSP) 1 



(51) 



x 



1 



in which x($ SP ) and Wf p ($ SP ) are calculated from (|4 
for a given fraction <j> SP . 



Calculation of the self-energy for a solution of Np 
noninteracting polymer chains 



These equations can be solved numerically and x(<f> SP ) 



SP 



< 



and Wi ($ 6 ) calculated for any values of < $ 
l,thus leading to the required W SP (k). As <& SP — ► 0, 
Wf p decreases to zero and also, from the structure of 
Eq. d48|), one expects a divergent behavior of Wf p as 
$ SP increases beyond a certain limit. The effective vis- 
cosity ?ysp of the medium containing the spheres and the 
translational friction coefficient Ct of the moving sphere 
are obtained from : 



?7sp - Vo 



Wf p ($ SP ) 
1 



Ct = Ct 1 = W SP (fc = 0) = 6nari (49) 



(1 + Wf P ($ SP )) 1 



Ii/2(^(1> SP ))K 1/2 (x(^P)) 



We consider Np Gaussian, noninteracting, free mov- 
ing polymer chains, each of length L = nl, dispersed in 
an incompressible, newtonian fluid of viscosity 770- Some 
external force F(r) generates the velocity flow. 

Once again, we will use the notation W for the self- 
energy of the fluid when only the polymers are present. 



Following the derivation in Sec. [IB, we can write the 



Navier-Stokes equation for the velocity field v(r) similar 
to (27) — but no spheres present — with the stick bound- 
ary conditions (|8|). Averaging the N-S equation over the 
configuration of the polymer chains and introducing the 
self-energy of the fluid W as in (Q), we can write the 
solution for the averaged velocity u(r) 
Eqs. (H, 0|l|)), where the force propagator Q 
given by : 



Q POL * F (see 

,POL • 



One can derive the above expression for the friction 
coefficient by applying the general analysis from Ap- 
pendix A to our particular case of a suspension of N — 1 
fixed and one uniformly moving spheres . Then, the re- 
lation between the total average force exerted upon the 
mobile sphere by the fluid and its velocity will be given 
by Eq. ( |A6| ) without the polymer term : 

Jdna(Q)^ = I^JdVLdQ! T(Q,,fl')\ u= (50) 

J dn rffi' Kr x (fi, o'A • u = Ct • u 



,POL 



>POL 



1 - kk 



?7ofc 2 



W POL 



(k) 

i-kk 



770(1 + w POL ){k 2 + Cj 



(52) 



POL 



^POL 



1 + W{ 



POL 



where, in the limit of k — > 0, the following approximation 
was made for the transverse part of the self-energy : 



W£ UL (k) 



-mtpoL 



ijo^r^Jic 2 (53) 



Note that now /C _1 is the generalized inverse of <?, to 
account for the presence of the fixed spheres in the solu- 
tion. Comparing then to the sphere contri bution to the 
self-energy from Eq. ([43]) and making use of (B7,B18) one 



can recover the formula for £ f displayed in ( [49 ) . When 
the filling fraction of the background spheres approaches 
zero, Wf p — > and one recovers the Stokes result for 
the friction coefficient of one sphere. 



Our goal is to find W POL (k) self-consistently, for any 
volume fraction <J) POi . If W POi (k) were the self-energy 
of the fluid containing Np — 1 chains, the contribution of 
one more chain added would be (similar to (143), but with 



„pol _ 



= 1/V and no spheres) 



W 



POL 1 



Np - 1 



n 

-Y. (expfik-^-S,-)]^. 



»>j>i 



8 



I' J 



(54) 



Here, in contrast to (f43|), 7C 1 signifies the generalized in- 
verse of the modified Oseen tensor Q because we imagined 
the last chain immersed in the effective medium created 
by the initial pure fluid plus Np — f polymer chains and 
the force propagator should be changed correspondingly. 

The calculations are rather laborious and are detailed 
in Appendix C. In the hydrodynamic limit of long wave- 
lengths, the expression of W POi (k) is : 



W POL (k) = --7 ?0 fc 2 $ POi (l + Wf OL )Q([3)t 

IT 

Q(^) = Ilog(l + A N ) + _L + £ (55) 

P V V 71 "/ V 71 " 27r 



Symbols have the following meanings: $ POL = 
is the volume fraction of the polymers ; R g 



4ttR-' 



,POL 



¥ fa 



the radius of gyration ; /3(<£> POL ) = ^j^Atc\ with £' the 
screening length in the effective medium, given by (52). 
In obtaining the previous relation we used the Kirkwood- 
Riseman approximation that amounts to consider only 
the diagonal terms in the Fourier expansion of Q~ 1 {s, s') 
(s is the arclength along the polymer chain), as detailed 
in Appendix C. We believe this approximation is phys- 
ically justifiable because we are working in the limit of 
long chains, when the off-diagonal terms are small. 

To get the self-consistent equations for the unknowns 
/3($ POL ) and W POL ($ POi ), we take the transverse part 
of W POi (k) by dotting in (1 — kk) and then we equate 
it with the expansion (p3|) to arrive at : 



/3($^) = 

w POL^POL, _ 



1 - f * POL Q(0) 



(56) 



D 



POL/^POL 



( t KR 

qpol^pol^ 



(57) 



1 



l + W{ 



POL 



= 1 - 



18 



POL 



where we employed Q(0) 



and Q 



(9ir^/4) Vo R g 



is the Kirkwood-Riseman friction coefficient for the non- 
free-draining limit. 

Limiting values for D POL and ?7pol can be obtained 
for small and large polymer filling fractions <i> POi : 



<S> POL « 1 



?7POL - VP 

Vo 

D POL -> 1 



r oi (large) < 1 



Vpol - Vo 



24 , 1 
^ 9 V 



oc as § POL 



(58) 



0.309. 



D 



POL 



(59) 



We remark that for low <j> POi the result for the vis- 
cosity is twice that of Kirkwood-Riseman |L4[ for one 
chain because we neglect the rotational terms in calcu- 
lating W POL (k). It is also notable that the constraint 
of the positivity of W POL in this effective medium ap- 
proach leads to the divergence of the viscosity and of the 
friction coefficient, which occurs at Qpfy — 0.309. This 
divergence is attributed to the rigid-body dynamics as- 
sumed for the polymer. 



IV. RESULTS FOR SPECIFIC EXAMPLES 



In this section we apply the above theory to two types 
of mixed polymers-spheres suspensions : 

S/P system : one mobile sphere among other fixed, rigid 
ones, in a polymer solution. 



which gives W POL as an analytic function of the vol- 
ume fraction $ POL . Also it is interesting to remark that 
(3(= = implies that the hydrodynamic screening is 
absent in the limit of k — > (observation that is actually 
manifest from the structure of (|55|)). The function Q(/3) 
is defined in Eq. (pq). 



We can calculate the viscosity of the effective medium 
?7pol (using (|l9|)) and the diffusion coefficient D POL of 
a polymer chain as follows : 



?7pol($ p ^) 
Vo 



m = wf OL 



P/S system : free moving polymer chains in a suspen- 
sion of fixed spheres. 

In addition, to gradually develop the theory we start 
by explicitly investigating two particular cases of the gen- 
eral models just mentioned, namely : 

1S/P system : one mobile sphere inside a polymer solu- 
tion. 

1P/S system : one polymer chain moving in a suspen- 
sion of fixed spheres. 
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A. One sphere moving in a suspension of polymers 

In this example we consider a sphere ol radius a moving 
uniformly with the velocity Vo inside a suspension of Np 
non-interacting, Gaussian, free moving polymer chains of 
length L. The flow of the fluid is generated by the mov- 
ing sphere. The motion of the particles is coupled to the 
motion of the solvent by stick boundary conditions. We 
also mention that the total force and torque acting upon 
each polymer chain are zero. 

Similar to Eq. (p7j), the equation of motion for the fluid 
and the boundary conditions can be written as : 

NP,n 

- 7? Av(r) + Vp(r) = F (r) + ^ 6(r - TL al )<r m (60) 



translational friction coefficient of the moving sphere is 
just the Stokes result : 



Ct = 67rarf ff 1 



(66) 



with T] cS = ?7o(l + Wf°^}. The translational diffusion 



coefficient as defined in ([21]) reads then 



D(<S> POL ) = 



6irarjo 



= 1 



18 



<I> 



POL 



(67) 



It is noteworthy that in this case there is no size cou- 
pling (no t — dependence) between the dynamics 
of the polymer chains and the dynamics of the moving 
sphere. The diffusion coefficient vanishes and the viscos- 
ity diverges at $ POL ~ 0.31. 



V • v(r) = 

Rqi = v(R Qi ) 



(61) 



B. One polymer chain immersed in a suspension of 
fixed spheres 



in which Fo( r ) is the force exerted by the mobile sphere 
upon the fluid given by : 



(62) 



F (r) = / dfl 5(r - R )er(r!o) 



Averaging the Navier-Stokes equation over the distri- 
bution of the polymer chains and introducing the self- 
energy of the solution : 



' TV P.i 



£(r,r')=W POi (r,r') = ( £ 5(r - R ai )a ai ) (63) 



\a,i>l 



we find the formal solution (following Eqs. ([l2 18)) 



v(r) = £*F | r 

900 



rj k 2 



W POL 



(k) 



(64) 



and we succeeded in transforming the initial problem in 
the Stokes motion of a single sphere in a fluid of viscosity 
Vpol- Because only the polymers are contributing to the 
effective viscosity of the solution, rf s is given by (|l^,[57]) : 



„cff 



Vo 



VO = W PQL 



18 ^PQL 



(65) 



As it was shown in the previous section — Eq. (|5C 
the hydrodynamic screening is absent for a suspension of 
polymers. In this problem, the sphere does not contribute 
to the self-energy of the fluid, which implies that the hy- 
drodynamic interactions are still unscreened. Then the 



This time we consider a free polymer chain of length 
L, immersed in a suspension of N5 rigid, fixed, uniformly 
distributed spheres of radius a. There are no interactions 
other than hydrodynamic and the flow is generated by 
some external force field F(r). 

The velocity field is described by another variant of 
N-S Eq. (^?J), with only one chain present : 

n 

- ijoAv(r) + Vp(r) = F(r) + ^ <J(r - R 4 )<r 4 

i=l 
Ns r 

+ J2 dn b 5(r-R b )<T b (68) 
b=i 

where all the symbols have the known meanings. 

Taking the configurational average of the equa- 
tion above and introducing the self-energy tensors 
S POL (r,r') and W SP (r,r') (see @) related to the in- 
fluence of the polymer chain and of the spheres on the 
fluid, we can cast the Navier-Stokes equation for the av- 
eraged velocity u(r) = (v(r)) in the following form : 

- i)oAu(r) + (Vp(r)) - W SP * u| r = F(r) + S POi * 

(69) 

the formal solution of which being (see Section II) : 

u(r) = g*F\ r + g*x POL *g*F\ r + ... (70) 

Here * is the usual convolution operator and Q{r, r') 
is the effective Oseen tensor modified to account for the 
influence of the spheres. Its expression is given by ( |4"6| ) 
using the expansion for small k (ESI). In this way our 
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starting problem was reduced to studying the station- 
ary dynamics of a polymer chain in an effective solution 
where the fo rce p ropagator is Following the re- 

sults of Sec. W%— Eqs. @ |f, M ©—and noticing 
that we must make the changes c 



= i, G -> and 
/C _1 , we find the solution of S POi (k) in first 



order in the number of scattering events 



^POL 



k) = T ^(exp[ i k.(S,-S ) )])J(K, 1 ) 



(71) 



Note that A3 1 is the generalized inverse of Q 

n 
3=1 



(72) 



Remembering that the effective medium has replaced 
the spheres and the pure fluid, we can repeat the 
derivation from Appendix C to get an explicit form 
for £ POL (k), with the changes W POL -> Wf p and 



„pol _ 1 



iPOL 



(k) 



-Vo- 



(1 + Wf p )Q(/? SP )fc 2 1 (73) 



3 V 



Psi 



Rg 



SP 



where £' sp and Wi were defined in Sec. Ill A 



To evidentiate the important dimensionless variables 
of this model, we introduce the parameters t — —j*- and 
xsp = 7^— and rewrite X POL (k) as : 



q 3>sp 

S POL (k) = --no — +3h 1 nrSPfiiSPv, „ , , ,.2 

7T N5 



t*[l + W? p (3> bP )]Q(tx S p)k 2 l 

(74) 



Both x S p(& SP ) and Wf p ($ SP ) are characteristic to the 
problem of the spheres immersed in a fluid and can be 
calculated from the system (Eg). 



In a similar manner, the diffusion coefficient of the poly- 
mer chain in the effective medium of spheres can be ex- 
pressed as : 



D(<P 



SP\ 



30? 



1 



P(tasp) 



(POL^SP) 



4 1 + W( 



SP 



tx S p 



(75) 



in which we used (|C24|), with W POL -► Wf p and 
( 3 — > Psp- The function P(x) is defined in Appendix C 



The effective viscosities, total and relative (which mea- 
sures the increase in the viscosity due to the polymer 
chain) are calculated from : 



off 



off 



= W- 



SP 



lim 



1 d 



^POL 



k^o r) 



<% 2 ± 



(k) 



W- 



sp , 9$ SP +3/ 



t a (l + W r f p )Q(tx S p) (76a) 



?7sp 



7T N 



7T N 5 
1 d 

fc^o rysp 

t 3 Q(tx S p) 



= -l™— ^M^(k) 



9$ SP 



(76b) 



with T] SP = 770(1 + M / i SP ), and isp($ SP ), Wf p ($ 5P ) 
obtained from (pig), as already mentioned. 

As will be shown below, Eqs. (|7|^) are the limit- 
ing formu las of the general situation presented in Sec- 
tion IV D, where the dependence of D and rf on <E ,SP 



and t will be discussed also. 



However, some remarks can be made. The arrest of 
the chain motion occurs at $fjjy — 0.49 — as found in 
Rcf. [ fill for a sphere moving in a random array of fixed 
spheres. The curves seem physically plausible. When 
the radius of gyration R g of the polymer chain is smaller 
than or comparable to the radius of the sphere (t < 1) 
the viscosities (both total and relative) are relatively in- 
sensitive to variations in the t parameter. Still, although 
being related only to the presence of the polymer chain, 
the relative viscosity depends on <& SP and t variables 
through Q(t xsp) (see [76]), so the contribution of the 
chain to the viscosity of the solution is a function of the 
effective solvent. 



One mobile sphere in a suspension of fixed 
spheres and polymers 



We now investigate the stationary dynamics of a 
suspension containing N5 uniformly distributed, rigid 
spheres of radius a and Np Gaussian, free moving poly- 
mer chains of length L. All spheres but one that is mov- 
ing with the constant velocity Uo,thus creating the flow, 
are fixed, and no interactions beside the hydrodynamic 
coupling occur. We calculate the diffusion coefficient D 
of the moving sphere and the total rf and the relative 
-2 1 effective viscosities as functions of the volume 



fractions <5> SP for spheres, <& POL for polymers and of the 
coupling parameter t = —. 

We proceed closely to the derivations exposed in the 
previous sections. First, by properly identifying the 
C/(k) tensor, we construct an effective medium replac- 
ing the polymers and the pure fluid, in this way reducing 
the problem to that of a self-consistent computation of 
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the viscosity and diffusion coefficient of an ensem ble o f 
spheres, as discussed in Refs. Jl^,||,0 and in Sec. [II A . 
Solving the resulting equations, we obtain the desired 
quantities, namely the diffusion coefficient of the moving 
sphere and the effective viscosities. 

The velocity field obeys the N-S Eq. @, with F(r) -> 
Fo(r) now being the force exerted by the moving sphere 
upon the fluid. Averaging first over the distribution of 
the polymers and then over the position of the fixed 
spheres and introducing the self-energies W POi (r, r') 
and S SP (r,r') as defined in @, we deduce the N-S 
equation for the averaged velocity u(r) = (v(r)) : 

- i) Au(r) + (Vp(r)) - W POL * u - £ SP * u = F (r) 

(77) 



where £' is the total screening length in the resulting ef- 
fective medium that replaced the initial pure fluid, the 
polymers and the spheres. Note that the poly mers do 
not contribute to the screening, as found in Sec. 



IIIB 



It follows then that the total and the relative (due to 
the presence of spheres only) effective viscosities of the 
solution and the translational friction coefficient of the 
moving sphere (see Sec. [ill A ) can be computed from : 



V 



cir 



m 

efi 



1 = W( 



POL 



Ef p 



d 



Vpol 



lim — 

k^o 77POL ok 2 



^SP 



(k) 



1 + W{ 



POL ' 



1 



sp^l p (k = o) = m-sp^sp 



^SP 



(83) 



As usual, we denote the self-energy of the fluid result- 
ing from the presence of the polymers with the W POL 
symbol and not with S POL to stress that this is the 
background effective medium we have related our cal- 
culations to. This also means that W POi depends only 
upon the polymer properties, but E — the contribu- 
tion of the spheres to the self-energy of the fluid — could 
depend on $ POL and on some coupling parameter, like t 
and,implicitly,on other characteristics of the background 
effective fluid. For W PO L we will use the explicit expres- 
sions found in Sec. 



IIIB 



(Eqs. B5LB6 



The formal solution of the N-S equation is then 
u(r) = £*F| r 



(78) 



with the modified Oseen tensor having the following 
form : 



0(k) 



1-kk 



77o fc 2 - W POL (k)-S^ p (k) 



(79) 



Next step is to find S SP (k) self-consistently. This can 
be done exactly as shown in Appendix B,but taking care 
of the effective medium that replaces the polymers, by 
using the new Q, substituting Ei for Wf p and chang- 
ing the meaning of which will now represent the total 
screening length. Thus : 

£ SP (k) = -6^0^(1 + Ei) ' 



Ii/ 2 (£)Ki /2 (£) 
|r7 (l + Ex) $ SP fc 2 Z(^) 1 + (longit. part) zz (84) 



1 



8 



Ii/2(j/)Ki/ 2 (j/) 9I 3 /2(y)K 3/ 2(y) 5Ii /2 (y)K 5/2 (y) 



Comparing this expression with (|8^) and using the def- 
initions of Ei and £ from © and of W[ OL from @), 
we finally obtain the system of equations for S SP (k) : 



Working in the k — > limit, we approximate S SP (k) 
with : 



s SP (k) = (- Vo ^ P - m^k") i 



(80) 



where £sp is the screening length felt by the spheres in 
the effective medium that supplants for the polymers and 
the pure fluid. Substituting back in the Oseen tensor 
formula and replacing W POL with its expression from 

(P|,P|) we g et : 



<?(k) = 



1-kk 



1-kk 



jjdCL + EOOfc 2 + r) eS (k 2 + £' ) 



(81) 



in which : 



Ei = W POL + Ef p 
£/-2 = 4sp 



1 + Ei 



(82) 



Ef p ($ SP ,$ POi ) 



= -$ 



18 
SP 



1 



1 



$ SP Z(x 



Ii/ 2 (a;)K 1/2 (a;) 



- 1 



(85) 



with x = x(Q SP ) = £ and the function Z(x) being de- 
fined in the previous equation. One could remark that 
there is no polymer dependence in the total screening 
length 

Solving this system for x($ SP ) and Ef p (<& SP , $ POL ) 
allows us to compute the diffusion coefficient of the mov- 
ing sphere from fl83) : 



D ^sp^po L) = ^m 

x(l 



I 1/2 (x)K 1/2 (x) 



(86) 



^L<S>pol ){1+ Ksp z{x)) 

TT\/TT 2 
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and the total and relative effective viscosities as 
?7 off 1 



1 (87) 



■ir 



Vpol 



1 



1 



1 + W( 



POL 



SP 



1 



\<S> SP Z(x) 



In Figs. §,|; [|§; | the results for D, the total 
i] eB /r]o — 1 and the relative 77 cfI /r/poL — 1 viscosities have 
been plotted against $ POL (first plot of each pair and 
Fig. |^) and against & SP (in the second plot). 

We have already seen that the total screening length £' 
as given by ( [So]) is function only on <J> SP (as the polymers 
were absent). It is also notable that there is no depen- 
dence upon the coupling parameter t = —. Both obser- 
vations are related to the absence of screening in a disper- 
sion of non-interacting polymers (as found in Sec. Ill B| ) , 
which prevents the spheres to sense the existence of the 
polymer chains. 

Two important points, at ^f/jy ~ 0.49 and at ~ 
0.31, where the viscosity of the solution and the friction 
coefficient diverge, are evident in Figs. ||-^|. Note that 
^mv + ^luv < 1' even if the polymer chains and the 
spheres are penetrable, non-interacting objects. 

At low $ SP and $ POL — > one recovers the Einstein 
result for the viscosity of spheres in a pure fluid and the 
Stokes friction coefficient for a sphere (D($" SP ,0) — > 1). 

When $ 5P — > (recall <& SP = Ns y 1 ) 1 we retrieve the 
results from the last section for one mobile sphere im- 
mersed in a polymer solution (compare Eqs. ( |S6|]8"lj ) to 
Eqs. ( ]67j , |65| )). The results of this limit are marked with 
the + symbol in Figs. || and ||. 



solution, for any volume f ractio n $ POL , following the 
procedure detailed in Sec. IIIB . Conceptually, we are 
modeling this system with a superposition of two effec- 
tive media, a fixed background related to the influence 
of the spheres and a foreground related to the polymers, 
coupled through hydrodynamic interactions. 

Apparently, there should be a symmetry between the 
present system and the reversed one (spheres in a so- 
lution of polymers), but this assumption proves to be 
incorrect, as we will further show. The main difference 
is that in the present case the size coupling is active and 
t = controls both the mobility D of the chains and 
the viscosity of the dispersion. 



We start with the Navier-Stokes (N-S) equation (J27 
for the velocity field v(r) in the suspension, with F(r) 
some external force density. After we perform the config- 
urational average ( |32| ) and we replace the averages over 
the forces with the W SP (r,r') and E POi (r, r') as de- 
fined in ([mJ) the N-S equation for the average velocity 
field u(r) is : 



%Au(r) + (p(r)) - W 



sp . 



u- £ 



POL 



* U 



F(r) 



where, as discussed in the previous sections, we substi- 
tute W SP for E SP to distinguish between the back- 
ground effective fluid associated with the spheres and 
the contribution to the self-energy of the added polymer 
chains. W SP (r, r') tensor is function only o f $> SP volume 
fraction and it was calculated in Sec. Ill A, Eqs. ( [47| , fl8| ) . 

To obtain the self-consistent equations for E POi we 
proceed as in the last section. First, we write the formal 
solution of the N-S equation as a convolution : 



D. Np polymer chains moving in a suspension of Ns 
randomly distributed fixed spheres 



In this section we explore the reverse of the previ- 
ous model. There are Np non-interacting, free-moving 
polymer chains of length L dispersed in a suspension 
containing a random array of N5 rigid, penetrable, 
fixed spheres. We wish to compute — using the effective 
medium approach — the translational diffusion coefficient 
D for a polymer chain and the total and relative (the 
spheres contribution) effective viscosities as function of 
<j> POi 7 $ SP and of the coupling parameter t = R g /a. 
The velocity flow is externally generated (one could con- 
sider it to be produced by the motion of one polymer 
chain without any influence on the final results). 

Once again, we develop the necessary theory in two 
stages. First, we substitute the spheres and the pure 
fluid with the corresponding effective medium presented 
in Sec. Ill A, then we calculate self-consistenty the trans- 



port properties of the polymers dispersed in the effective 



u(r) = g * F 



(89) 



where the modified Oseen tensor has its Fourier trans- 
form given by : 



C?(k) = 



1 - kk 



770 fc 2 - VFf p (k)-E POi (k) 
1 - kk 
~ 77 (1 + E 1 )(fc2+£'- 2 ) 

Ei - Wf p + E POL ; £,'- 2 - isp + ip2 ° L 



1 + Ei 

in which we inser ted Wf p from © and we used the 
following approximation , in the k — > region : 



(k) = (-vo&oL ~ V^i° L k 2 ) 1 



(90) 



Also, £' is the total screening length and £sp is function 
only of the concentration of spheres. We can readily write 
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the solution for E POi (k) by observing the structures of 
the Oseen tensors found in ( |52[ ) and (|90|). In compari- 
son with t he effective medium theory for polymers (as in 
Sec. [II B and Appendix C), we need to transform (|55| ) 
by working with £' instead of £' SP and with the total 
self-energy Ei instead of W POL . We immediately arrive 
to : 



iPOL 



(k) 



-77o$ POL (l + I]i)Q(/3)fc 2 l (91) 

7T 



with = -|f and the function Q(x) defined by (55) 
Comparing to (|90"|), we obtain the system : 



^POL _ u ^POL , 



... _ _ $^ (i + + Ef uv ) Q(/3) (92) 



Introducing the natural dimcnsionless variables 



SP\ _ a 



through 



and t — we finally determine E ° (k) 



E p ($ , $ 6P ,i) = [1 + Wf P ($ 6P )] 

1 



(93) 



1 - % $ POL Q(/3) 
/3 2 = (ix) 2 (l--$ POL Q(/3)) 

where Wf p and x($ SP ) are computed from the system 
( [|| ) as functions of <S> SP only. Note that now /?(= R g /C) 
is in general a non-zero function of $ POL 5 $ 5 ' p and 
t, which means that the polymers are subjected to a 
screened hydrodynamic interaction, in contrast with the 
previous system. 

The effective viscosities — total and relative, due to 
added polymers — are determined from : 



off 



- 1 = 



W SP + E POL 



(1 + wf p ) 



1 - f <f POi Q(/3) 



Vsp 



1 = 



l + W{ 



1 

SP 



^POL 



1 



1 _ l$POLQ (/3) 



(94) 
- 1 



For a mobile polymer chain the translational diffusion 
coefficient (actually defined as relative to the Kirkwood- 
Riseman result) is shown to be (see Appendix C) : 



D 



KR 



with C/qq 1 being the zero-th coefficient of a Fourier series 
expansion of (Q~ ) (see Appendix C). Because in our 
approximations Q qq i is diagonal, Qqq becomes ( C23 ) : 



(Go) 



3TT7] Q R g 



P 



P(/3) 



(96) 



The diffusion coefficient of any moving polymer chain is 
given by : 



D 



3^ 1 PQ9) 9 POL 
A 1 _i_ uiSP a v 1 „ w 



4 1 + Wf p /3 



J Q(/3)) (97) 



The functions P (:r) and Q (x) have been defined in Ap- 
pendix C, Eqs. ( |C23| , |C28| ). 



The prominent feature of the final results for the so- 
lution viscosity and the diffusion coefficient of a polymer 
chain ( p^ , p7| ) is their dependence on the t parameter — 
the stationary dynamics of the spheres and the polymer 
chains is size coupled — as advertised in the beginning. 
Also note that in contrast to the problem of spheres in 
a polymer suspension, the variable f3 = R g /t;' depends 
upon $ SP too, which implies — Eq. ( |9~4|) — that the rela- 
tive contribution of the polymers to the viscosity is en- 
hanced, function of $ SP and t, by the presence of the 
background effective fluid related to the spheres. It is 
sensible to believe that these new phenomena have been 
generated by the spheres acting as fixed scattering cen- 
ters of the velocity field and had we analyzed a system 
of polymer chains and Brownian spheres, we might have 
observed that the dynamics of the two types of particles 
would have been decoupled, in the mentioned sense (no 
t dependence, /3 function only on <J> POi ). 

Checking the limit $ SP — > we recover the results 
for the suspension of polymers only (Wf p — > 0, /3 — > 0, 
compare to (p7|)), and in the limit $ POL — > we re- 
trieve the results for the 1P/S problem : j3 — > txsp', 
<f> POL -> x747ri?3/3 = ^-t 3 (compare to Eqs. @ |76|)). 

The values of the total and relative effective viscosities 
and of the diffusion coefficient given by Eqs. have 
been plotted as functions of $ POi and <& SP at different 
t. Four sets of data corresponding to t — {0.01; 0.1; 1; 10} 
have been calculated and displayed. 

In Figs. Jt|-|io|1 the diffusion coefficient D of a mov- 
ing polymer chain is plotted against $ SP with $ POL pa- 
rameter and for each of the t values mentioned, and in 
Figs. |l2|-|l4| D is plotted as a function of $ POL with § SP 
parameter, with t — {0.1, 1, 10} for convenience. 

D(<fr POL ) depends linearly on $ POL and for fixed 
$ SP and $ POL ,the diffusion coefficient D decreases as 
t increases — longer chains are less diffusive — so the t pa- 
rameter acts like a mobility selector. The sensitivity of D 
to t variations is higher when the size (R g ) of the chains 
(95) becomes comparable to the radius of the spheres (com- 
pare the relative drop in D at the same $ SP in Fig. 12 
(t = 0.1) to Fig. |l| (t = 1), and Fig. |l| to Fig. [L4 
(t = 10). 

The same observation can be made when looking at 
the <J> SP dependence of D, in Figs. 0-0, for the same 
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values of t. Additionally, for t — {0.01; 0.1}, D exhibits a 
change in the curvature and an almost plateau region in 
the middle of the $ SP range for higher <& POL (Figs. ^,||). 

We now consider the behavior of the limiting case cor- 
responding to the 1P/S system as shown by the + sym- 
bols in Figs, f? -10 1- One can notice a change in the cur- 
vature of the diffusion coefficient of one polymer chain 
when going from t = 0.01 to t = 10. It appears that 
for relatively small (in comparison to the spheres) poly- 
mer chains (t = 0.01) or relatively large (t = 10) a 
plateau region appears, where D is slowly varying with 
respect to $ SP . Physically, this could signify that when 
the polymer is small (large), a certain built up in $> SP 
is required before the diffusion coefficient D is affected 
more markedly. Additionally, for a given radius a of the 
spheres, the difference in the mobilities of two polymer 

chains of distinct lengths (ii 7^ ti\ti = — s —) is maxi- 
mized for a certain range of $ SP values, fact depicted in 
Fig. [fl] , where the difference in the diffusion coefficient 
of one polymer chain for four pairs of t values is plot- 
ted against the filling fraction $ SP (e.g. $ SP € [.15; .3] 
roughly for two types of polymer chains with t\ = 0.01 
and t 2 = 1). This may have applications in the design of 
polymer sieves. 

Comparing now the behavior of the polymer moving 
in the suspension of spheres (Figs. [0-|O[) to that of 
the mobile sphere inside a polymer dispersion (Fig. |^) 
at equivalent relative sizes of the tracer particle (t = 10 
in the first case is similar to t = 0.1 in the second one), 
we notice that the D($ SP ) for a small mobile polymer 
chain (t — 0.01; 0.1) is higher than that of a small sphere 
(t = 10) moving in a polymer suspension, but a larger 
sphere (f = 0.1) exhibits an increased mobility (D) com- 
parative to a large polymer chain (t = 1; 10). Eventually, 
because ($f> p v ~ 0.49) > — °- 31 > tne polymer 

immersed among spheres will retain a non-zero mobility 
at higher volume fractions than the sphere moving in a 
polymer suspension. 

The same discussion applies in the general P/S case. 
The mobility of a small t = 0.01 or large t = 10 polymer 
chain is only slightly influenced by the volume fraction 
of the spheres (see slopes of graphs in Figs, [ijjn] in the 
middle region). This feature indicates that, for a given 
radius a of the background spheres and a pair of polymer 
chains characterized by distinct t values (e.g. t\ = 0.01 
and £2 = 1), the difference in the diffusion coefficients of 
the two chains will reach its maximum for a certain range 
of the spheres volume fraction <j> SP . 

The two points of divergence — 0.31 and 

$f, P y ~ 0.49 are clearly marked. Because $£7^ < $f> P y 
the suppressing of the mobility of the moving particles 
(either sphere or a mobile polymer chain) is stronger 
when polymers are added than when the sphere concen- 
tration is increased. 

Comparing the data for D($ SP ) in Fig. ||, for the S/P 



example (note that there is no t dependence, to those for 
D($ POL ) at t = 0.1 in Fig. [if and t = 10 in Fig. ||, the 
same observation as in the 1P/S case can be made : at 
the same t ratio between the sizes of the moving object 
and of the background component and for the same back- 
ground concentration, a small (t = 0.1) polymer chain in 
the P/S system is more mobile than a small (t = 10) 
sphere moving in the S/P system, but as t gets larger, 
a big sphere becomes more diffusive than a long polymer 
chain (t = 10). In this situation, the control parame- 
ter is the concentration of the added (e.g., in the P/S 
example, polymers are the added elements) components, 
but the same conclusion is reached if comparing the data 
for _D($ POL ) from Fig. H in the S/P case to those from 
Fig. g (t = 0.1) and Fig.]l0| (t = 10) for the P/S system, 
the control parameter now being the concentration of the 
background elements. 

The limiting case of one polymer chain (Np — > 1) 
is very well approached in the P/S example, as one 
can notice in Figs. f7]-[l0| by comparing the uppermost 
curves (1P/S case) to those corresponding to the pa- 
rameter value $ POL = 0.01, representing D($ SP ) when 
t = {0.01; 0.1; 1; 10}, respectively. 

Viscosities calculated from (^4[) have been plotted in 
Figs. [l5|-|2ll : the total effective viscosity change in 
Figs. |15|jl6| as a function of $ POi and against $ SP in 
Figs. [l^Jlq with t — {0.01; 10}; the relative effective vis- 
cosity in Fig. [l9] as a function of $ POi for t = {1; 10} and 
as a function of $ SP in Figs. [|], |2l] with t = {0.01; 1; 10}. 
Two observations could be made regarding the viscosi- 
ties : 

a) both for 77 cff (<I> POi ) with $ SP parameter 
(Figs. |l5|,|T6|) and for r] cS ($ SP ) with $ POL parameter 
(Figs. " , p^ j, as t increases, the total viscosity change 
is more sensitive to variations in $ POL or $' SP (com- 
pare Fig. [l5] when t = 0.01 to Fig. [l6| when t = 10 and 
Fig. (* = °- 01 ) to Fi §- HU = 10)). We note that a 
larger polymer chain ( in comparison to the radius of the 
background sphere) is more sensitive to the environment. 
This behavior clearly contrasts that encountered for the 
total effective viscosity in the S/P case (see Fig. 

b) the relative effective viscosity due to the presence 
of the polymers (Figs. [l9| fii"|) depends not only on $ POL 
but on $ SP and on t as well, once again in contrast to the 
S/P system (see Fig. ||). As suggested before, t acts like 
a control parameter that when increased (compare t = 1 
to t = 10 calculations in Fig. or Fig. |20] to Fig. |l]) 
makes the contribution of the polymers to the viscos- 
ity more markedly dependent on both <f> POi and $ SP . 
Actually, when the radius of gyration R g is comparable 
to the sphere diameter a it < 1) , the relative viscosity 

1 is relatively insensitive to variations in <& SP or t 

'JSP f —, J . . 

(compare Fig. |20| to t = 1 data in Fig. |21| ). 

Obviously, these observations may have some signifi- 
cance from an experimental point of view. 
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V. CONCLUSIONS 

In the present paper we have investigated the steady- 
state transport properties of a heterogeneous suspension 
of rigid polymers and spheres without excluded volume 
interactions, applying the multiple scattering formalism 
[HHID and the effective medium theory @,|5|Jl|. 

The main approximations and features of our approach 
were : 

a. the absence of any interactions among polymers 
and spheres. All particles are penetrable. 

b. the solvent is assumed incompressible and it is de- 
scribed by a linearized Navier-Stokes equation. 

c. the stationary, time-independent limit is assumed. 
All results are valid in the long range hydrodynamic 
limit (k — > 0). 

d. the dynamics of the solvent, polymers and of the 
spheres are coupled by stick boundary conditions. 

e. the spheres are monodisperse, having radius a. 

f. the polymers are monodisperse with chain length 
L, and they translate freely as rigid bodies. We 
do not consider the rotation of the polymer as a 
rigid body. For convenience, the distribution func- 
tion for the inter-segment distance of any chain was 
taken to be Gaussian. 

g. S(k), the self-energy of the fluid in the presence of 
the immersed spheres and polymer chains was ap- 
proximated with the first-order component E' 1 ' (k), 
in an expansion in the number of distinct scatter- 
ing events (one event means one interaction of the 
velocity field with one particle). When no corre- 
lations (interactions) among particles occur, this is 
actually the leading term of the exact S(k). 

Under these approximations, two general systems have 
been analyzed : 

1. Moving sphere in a solution of fixed spheres and 
freely translating polymer chains. 

2. Moving polymer chains dispersed in a suspension 
of fixed spheres. 

Also, two significant limiting cases of the previous models 
have been investigated : 

a. Moving sphere in a solution of polymers. 

b. Moving polymer chain immersed in an array of 
fixed spheres. 



For each example the underlying theory has been exposed 
and developed and equations were explicitly displayed, 
from which the translational diffusion coefficient D of the 
probe particle and the effective viscosities of the solution, 
total rj cS /rjo — 1 and relative, have been computed. 

We have identified the proper dimensionless variables 
characterizing the problem as follows : j3 = -gf , x = ^ 

and t = — where R g and a are the radius of gyration 
of a polymer chain and the radius of a sphere, and £' is 
the total screening length in the solution. 

It has been shown that the polymers do not contribute 
to the screening of the hydrodynamic interactions at 
large length scales and, as functions of polymer filling 
fraction $ POL , both the effective viscosity rf s and the 
translational friction coefficient £ t diverge, in all exam- 
ples, at a volume fraction $div — 0-31 much less than 
the corresponding divergence point for the $ SP depen- 
dence (<f>ls P IV ~ 0.49). While for suspensions of uni- 
form particles showing no long-range order — in partic- 
ular spheres — there is wide experimental support for the 
divergence of viscosity and of the friction coefficient for 
®d P iv e [0-5,0.6] (e.g. Ref. f|), in the case of solu- 
tions of rigid polymers and specifically for suspensions 
of polymers and spheres less data are available, to our 
knowledge. 

Perhaps the most interesting parameter is t, which is 
controlling the response of ry eff and D, both for spheres 
and polymers, to changes in the environment (variations 
in $ SP and Q POL ). The understanding of the role of the 
t variable might be of interest in experimental applica- 
tions, mainly related to gel permeation chromatography 
and macromolecular sieving. A rather noticeable predic- 
tion of our analysis related to the process of diffusion of 
the polymers in a suspension of fixed spheres is the ex- 
istence of an optimum range of the volume fraction $ SP 
of the spheres that maximizes the difference in the diffu- 
sion coefficients for two ty pes of chains characterized by 
distinct t values (see Fig. [n] for the 1P/S system and 
Figs . 0-0 in the P/ S case) . 

It has been also noticed that at the same t ratio be- 
tween the size of the tracer (moving) particle and the 
size of the background clement and the same background 
concentration, the spheres and the polymers are behav- 
ing differently. When the polymer is small in compar- 
ison to the spheres, it is more diffusive than a small 
sphere moving among polymer chains, but a large poly- 
mer chain in a suspension of fixed spheres is less mobile 
than a large sphere moving in a polymer dispersion. We 
must note that this observation is valid before reaching 
^-0.31. 

Despite the approximations made to yield the calcu- 
lations analytically tractable, the formalism presented is 
essentially general and also applicable when excludcd- 
volume (or other types of) interactions are turned on, 
although in practice the effort could be quite laborious. 
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Yet, in this straightforward approach we have managed 
to draw a physically reasonable picture of the dynami- 
cal behavior of such heterogeneous solutions of polymers 
and spheres and some conceivable findings emerged. Fur- 
ther extensions of this work are possible, for example the 
study of rigid rodlike polymers and spheres, microemul- 
sions (deformable particles) or investigating polydisperse 
suspensions of spheres and polymers. This last aspect 
is particularly meaningful in industrial applications (e.g. 
ceramic fabrication, where processing at volume fractions 
higher than the limiting values Qjjjy, &r>iv * s necessary 
but very difficult due to the divergent viscosity), because 
experimentally it has been shown that a proper adjust- 
ment of the particle size distribution could induce sub- 
stantial increases in the limiting volume fractions |l9| , 
but theoretically this observations is still not fully un- 
derstood. 

We seek to address this question and the physically 
more appealing problem of a heterogeneous suspension 
of particles with potential interactions in a future paper. 
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the relative distance : G(r — r'). Then, Eq. 
rewritten as : 



v(r) = G * F + G r 



Gr : f, • (Tb 



can be 



(Al) 



The sums over i,j,... run from 1 to n,the sums over 
a, f3, . . . from 1 to Npand the sums over b, c, . . . types of 
indices run from 1 to N5. Inserting the previous expres- 
sion in the boundary conditions (Eq) we get : 



U a + U> X S ak = G * F| Q + Gakjij ■ Cp] + G a k,b ■ &b 







G*F| 



G, 



(A2) 



Here u a and u are the velocity of the center of mass 
and the angular velocity of chain a. To eliminate the 
unknown quantities cr a i and (Tb we need the generalized 
inverse operators K^ 1 and K^ 1 Jl^,[l5| defined by : 



^K- 1 (S,,S J )-G(S ctJ -S c 



.7 = 1 



J dQ' b ' K b \r, r") ■ G(r" - r') = 1 6(Q b - Q' b ) (A3) 

where r(r2b),r"(SlJ ) ') and r(Q' b ) are position vectors (about 
the center of mass frame) of distinct points on the surface 
of the sphere b and S a j is the position vector of the i-th 
bead of the chain a with respect to its center of mass. 
Applying the inverse operators in (A2) we obtain : 



APPENDIX A: SOLUTION OF THE 
NAVIER-STOKES EQUATION FOR A 
SUSPENSION OF N P POLYMER CHAINS AND 
N s SPHERES 



In Sec. [I B we have found that a formal solution of the 
Navier-Stokes equation (|2^) is given by Eq. (|3(]), where 
the force densities <r Qi and (Tb exerted upon the fluid by 
the j-th bead of the polymer chain a and by the point R5 
on the surface of the sphere b, are unknowns. Using the 
stick boundary conditions (|28|), we eliminate the forces 
to express the velocity field in terms of the single-object 
flow propagators T Q of any chain a and T& of any sphere. 

To simplify the notation, we replace the variables 
Rb,R Q ; with the indices b and ai, respectively, when ap- 
pearing in the arguments of the Oseen tensor G and we 
apply the Einstein convention of summation over repeat- 
ing indices. In the case of the spheres this will also imply 
a surface integral over the solid angle dQ. We drop the 
chain index a for relations generally valid independent 
of the chain variables. The convolution operation is de- 
noted by * — f * g = fdr'f(r — r') ■ g(r') with f and g 
tensors or vectors — and it is understood that its result is 
a function of some position vector r, when not otherwise 
specified by a subscript. Due to the translational invari- 
ance in the suspension, G(r, r') is actually a function of 



K <*Lk ■ (ua - G * F\ ak - G ak . b ■ <T b 



fc=i 



G akM ■ Vf3j\ 07 L a ) + K a Lk ■ ( W « X S "fe) ( A4a ) 



k=l 



(T b (flb) 



K 



/.// 



K fc V • G fc 



G*F| 



K 



bb' 



(A4b) 



Note that there is no summation over a in the first equa- 
tion and over b in the second one. 

Analyzing the structure of these force factors, it ap- 
pears manifest that the force upon the fluid due to one 
object, is proportional to the relative increase in the ve- 
locity of the fluid about the local velocity owing to the 
action of all other bodies and external forces, which is a 
sensible physical expectation. Also, the inverse operators 
. and K^} will produce the friction coefficients of 
the polymer chains and of the mobil e sp here, after one 
performs the configurational average (|32]) and sums over 
the indices i, k in the polymer case or integrates over Qb 
and fl' b for the sphere. 

As shown in Ref. |12] ], the ang ular velocity of the 
chain uj a in the force equations (A4) leads, eventually, to 
the appearance of the rotational and cross translational- 
rotational/rotational-translational friction coefficients of 
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the polymer chains. The former has been calculated in 
Rcf. |12| and the latter vanish upon preaveraging. In the 
present work we will ignore the rotational term, which 
does not alter the screening properties of a polymer so- 
lution. Then, using the conditions ( p9] ) that the total 
external force and torque acting on any polymer chain 
are zero, we get the velocity u Q : 



gr 1 • E K ^ • G * F U + gr 1 • E K ^U 

i,k i,k 
n 

■G a k,b ■ (T b + gt 1 ■ E K alak ■ G ak,l3] ' fftj 
i,k 

n 

g^E^ 1 : gt-gr 1 = i 



0^a 

(A5) 



Inserting back in (A4) and relabeling the dummy indices, 
we finally calculate : 

C ai — Tai,ak ' G * F| afe + T a i,ak ' G a k,f3j ' & f)j \ 0^t a 

i -L ai,ak ■ G a k,b ■ <Tb (A6) 
O'b = Tbb' ■ G * F\ b , + Tbb' ■ Gb',ai ■ Cai 
+ Tbb' ■ Gb'cCc\b^ c 

where there is no summation over a in the first equa- 
tion and over b in the second and we have introduced the 
notation : 



- ai,ak 



ai,ak / j ai,al of ctV ,ak 



U'=l 



T bb ' - -K b6 } 



(A7) 



Following indefinite iterations in (|A6[) and substituting 
the results for cr ai and <Tf, back in ( |Al[ ), we obtain the 
velocity field as : 

v(r) = G * F + G ra i ■ T ahak ■ G * F\ ak 
+ G rb ■ Tbb 1 ■ G * F| b , 

+ G ra j • T a i ia k ■ G a k,/3j ' Tpj^pilp^ ■ G * Fj^ 



+ G r b • Tbb' • Gb'c • T cc i ■ G * F 

+ G r 

+ G r b ' Tbb' ' Gb'. 



b^c 



(A8) 



16' 
I ak 



It is worth noticing that both operators still depend 
upon the centers of mass through R Qi = R° + S a i and 

Eventually, the above definitions allow us to write ( |A8|) 
in a physically more suggestive manner, as a sequence of 
multiple convolutions that are representing the scattering 
events : 

v(r)= G * F + E G * {T Q *G*F + E G *Tb*G*F 

a b 

+ EE G * { Ta *G*Tb + T 6 *G* T Q } * G * F 

(A10) 



a b 



+ EE G*T a *G*T /3 *G*F 
+ EE G*T b *G*T c *G*F + 

b c&b) 



where the dependence on r of the multiple convolutions 
is implied and the sum is continued over all permissible 
scattering sequences. 

This is Eq. (pfi) we wished to derive. 



APPENDIX B: SELF-CONSISTENT 
CALCULATION OF THE W(k) TENSOR FOR N s 
SPHERES IMMERSED IN A FLUID 



The main tool we need is the following expansion in 
spherical harmonics : 

oo I V 

A(r(n),r'(Q'))= E E E ^m;*wlW<W m ,(fi') 

l,l f >0 ra— — l m' — — l' 



A 



Next, the single chain and single sphere flow propagators 
T a and Tb will be defined as : 

n 

T Q (r - r') = E S ( r ~ R ») 5{r' - R ak ) (A9a) 



i,k=l 



T h (r-r' 



') = JJdQ b dSl' b S(r - Rb) Tbb' 5(r' - R 6 ,) (A9b) 



where r and r' are the position vectors for points on the 
surface of a sphere (|r| = |r'| = a) with angular direc- 
tions f2 and f2', Yi m (Q) are the spherical harmonics and 
A is any tensor (vector)-like quantity. 

Consider now Np — 1 fixed, non-interacting spheres of 
radius a immersed in a pure fluid of viscosity 770- One 
can calculate the exact self-energy S SP (k) of the sus- 
pension as in Sec. |ll B| (dropping the factors related to 
the polymer chains) but a more direct approach is to 
seek a self-consistent approximation , W££(k) , of the 
exact self-energy. As discussed in Sec. II A, we can re- 
place our system of spheres and pure fluid averaged over 
the distribution of the spheres, with an effective medium 
in which the force propagator is Q SP (k) given by Q). 
Adding one more sphere, its contribution to the approx- 
imate self-energy of the fluid is found from the second 
factor in ( p3| ) with the replacements c POL — > 1/V and 
K 1 xr^: 
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W ap (k 
Nc-1 



x exp[ik- (r(0) - r' («'))] (B2) 



where /C 1 is the single sphere generalized inverse of the 
modified Oseen tensor Q : 



dndn" /c _1 (ji, o") • g(r(n"),r'(n')) = t6(n-n') 



(B3) 



Expanding in (B2) in spherical harmonics we get : 

lm I'm' 

x exp[ik-(r-rO]Y Im (fi)Y?, m ,(fi') (B4) 

But the quantities of interest (shear viscosity, diffusion 
coefficient) are related to the transverse part of W SP (k), 
which should be isotropic in the k space (in our linearized 
model) . Thus we can perform the angular integral in the 
equation above by choosing k || z direction and employ- 
ing the plane- wave expansion : 

oo 

exp[ik(z - z')] = J2 ( 2/ + W + *) 

l,l'=0 

x ji{ka)3v{ka)Pi{cosO)P v {cos6') (B5) 
cos# = k-f; cos0' = k-r'. 

we arrive at : 

oo oo 

W SP (k) = -4ttc sp ^^(-1)^ w '^/(2Z+T) 



i=02'=0 



x ^{2V + l)j l (ka)j l ,(ka)Ki Q ) m (B6) 

Note that this equation is valid for any k values. A partic- 
ular value of the self-energy W S P (k) that is meaningful 
and readily obtainable from (|B4[), is : 



W SP (k = 0) = -^c SP K^. 



1)0 



(B7) 



For evaluating 1C l7 ^. llm ,, we start by expanding the def- 



inition (B3) in spherical harmonics to transform it in the 
equivalent form : 

oo I" 

^ ^ ^Tm;l"m" • &l"m";l'm' = A fa' 8 mm > (B8) 
Z"=0 m"=-l" 



Substituting £?(r,r') — Q(r — r') with its Fourier trans- 
form (calculated as in (||)) yields : 



Hind' m' 



]i(ka)] v (ka) 



-{~i) i i i+v i dkk 2 ^r^rrr^ p 



(BIO) 



with j/(Zca) the spherical Bessel functions of order Z and 
Fim:/'m' the tensor : 



P im; i'm' = / dfifc (1 - kk) YJ m (fife) Y,/ ro , (fife) (Bll) 



Now inserting VFf p (fc) from ( E^ ) and integrating over 
the k variable produces pa] : 



^4(5)^(5) 



" %a(i + wf p ) v +r? 

lm;l'm' 

(B12) 

Here, I and K are the modified Bessel functions of the 
first kind and £' is the effective screening length in the 
suspension of spheres from (^6|). 

The construction of the inverse K^m-i'm' rec l u i res an in- 
volved procedure , described in Rcfs. QQ. It is based on 
the observation that £7 — in its matrix representation — 
can be decomposed into a block diagonal (Z = Z') £/£> 
and an off-diagonal Qod part (Z ^ Z'). We could write 
then formally Q = Q d + Qod, which leads to, by tak- 
ing the inner product from left and right with 7C _1 and 
K.^ 1 , respectively, and iterating indefinitely : 



is~ — 1 — 1 is* — 1 /~> — 1 



(B13) 



We compute then first (/C^ 1 



>lm:lr. 



(K^) lm . !lm , = m a (1 + Wf p ) - — 1 — P^ 



(B14) 



in which the last factor is the generalized inverse of 
Pim ; /m' defined by the relation : 



where the coefficients of the Oseen tensor expansion in 
spherical harmonics are : 



JJdndn' e[r(n),r'(fi')] Yj m (n)Y Fm ,(Q') 



(B9) 



n" — — l 



(B15) 



and then one can calculate the non-diagonal part (Z 
«'+2)of/C^ m , : 
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lm;l' m' T / a \ t/ / a "\ lm;l'm' 



(B16) 



Finally, we return to (B6) to limit the expansion to fac- 
tors of order A: 2 and using also the constraint Z— Z' = {0, 2} 
we get : 

W SP (k « 1) ~ -47rc SP [K m ) 00 f (ka) + SK w ) w il(ka) 
- 2^5 j 2 (fca) j (ka) Re/C^ool (B17) 



The actual values of the operators ICi^-i'm' are : 

1 (B18a) 



_ 3 1 + Wf p 

^00:00 — o fl7 ?° 



2 Ii/ 2 (f)K 1/2 (^) 



K 



1 

10-10 — n a7 ?° 



1 + Wi 



9 '"I 3/2 («)K 3/2 (^) 



1 TT-l 



a \ ~W , I a \ 10;10 



(B18b) 



^20:00 - „ /7T7T a7 ?0 



1 + Wf p _ a 
^ — rarTTT — ttt ^oioo (B18c) 



30 '~ii/ 2 (f)K 5/2 (f 

with the 1L _1 tensors given by : 
JL^. 00 = 121 - 2zz 

IL 2 "o;00 = sJ^l + uJ^tt 



(B19) 



The last step in obtaining the self-energy of the sus- 
pension of fixed spheres is to expand the spherical Bessel 
functions in ( B17[ ) and, collecting all factors, we get : 



SP 



W* p (k) = -ftr^ W prr l 



3 
2 

2 

27 I 3/2 (f )K 3/2 (f ) 
1 1 



<f SP (1 + Wf p ) fc 2 
1 



Ii/ 2 (f)K 1/2 (f) 



— (121 -2zz) 



(B20) 



10V6 hMf )K 5/2 (f ) 



-1 + 24^-zz 



that reduces immediately to (E7|), the equation we set to 
derive. 



APPENDIX C: SELF-CONSISTENT 
CALCULATION OF W POL (k) FOR A 
SUSPENSION OF N P MOBILE POLYMER 
CHAINS 



a) change the variable to R(s) where s is the posi- 
tion of the z-th bead measured as the length of the 
arc along the chain, that is i = f, Z being the Kuhn 
length. 

b) convert the sum 53i=i("") the integral 
j Jyds (• • -),and the Kronecker delta &y to I6(s — 
s>). 

c) the Fourier transform of a quantity (A(s, s')) — the 
average being taken over the distribution of the 
segments about the center of mass of one polymer 
chain (assumed to be Gaussian) : 

(A(.y)>= £ A g9 ,exp(^-^ 

q,q'— — 00 ^ 

^ qq ' = T2 I ds [ ds'{A(s,s')) (CI) 
L Jo Jo 



x exp 



2iirqs 2iirq's' 



L 



L 



The generalized single chain inverse JC-^Sj.Sj) of the 
modified Oseen tensor Q a (Si, Sj) of the polymer solution 
is defined by the relation : 

n 

yj/C~ 1 (S,, Sj) ■ Q(S a j, S a k) = S ik l (C2) 
i=i 

Averaging over segments distribution and applying the 
previous transformations we get : 

^ f ds' K-^s, s') ■ g(s', s") = 6(s - s") 1 (C3) 
' Jo 



the Fourier transform of which being : 
+00 ^ 

£ K tf ■ g i'i" = ^ V' 1 

q'= — 00 



(C4) 



-1 _ 1 n -\ 



or, in operator notation JC 

Our starting point is Eq. (^ij) converted to an integral 
form : 

W POL (k) = -c POL i ( L ds \ L ds 1 (exp[ik-(R(«) 
' Jo Jo 

-Ris'Wii^^s')) (C5) 

- J2 is* 1 ) ■ [ [ dpdpl ( /C " 1 ( s ' p )) • <K"Vy)>} 



In order to compute the self-energy of a solution of 
polymers, we need to make the transformations : 



where the averages are taken over the distribution of the 
segments about the center of mass of the polymer chain. 
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For a Gaussian probability distribution function we have 
the following result jlo| : 



(exp[ik • (R(s) - R(s'))]) = exp 



6 



(C6) 



Inserting the above expression in (C5), Fourier expanding 
the (7C _1 ) factors and using their definition gives : 

W POi (k) = W POL (k) + W POL (k) 

w^(k) = - c ^ £ KAk)s; q H (era) 

q,q' — — oq 

1 Jo Jo 



q,q' u,u)' 

2i7TUJp 2i7ru f p f \ i 1 

K 

qu) '^uj'q' 



i*,(k) exp + niruL ] a 



with the notation 



1 /" L 








k 2 l\s-s'\~ 




/ ds ds' 
Jo 


exp 




6 




2inqs 


2inq' 


s'~ 




x exp 


L 


L 







(C8) 



In the previous relations we implicitly assumed that all 
tensors are multiples of the unit tensor, which will be 
shown to be true under certain approximations. 

The integrals over p,p' in the W POi expression are 
straightforward, yielding : 

Wr L (k)=c POL (^E^'W^o 1 ^! (09) 

qq' 



Now I*, (k) can be evaluated to be 



l *qq'i k ) ~ J2 



k 2 R 2 g 



2irq 





[m 


2 

- QQ 


1 






Ms) 2 ] 



(CIO) 



[l-CX P (-fc 2 i?2)]} 



in which R g = y Pjr- is the radius of gyration of any 

polymer chain. What remains to be calculated is the 
Q^ 1 operator. 

Approximating the transverse component of W POi (k) 
with (^3|) and using the translational invariance of the 
suspension we compute first the effective Oseen tensor : 



g(K(s),K(s'))=g(R s -R' s ) 

= f 0(k) exp Hk • (R s - R'J] 

Jk 

1-kk 



6(k) 



m k 2 - W POL (k) 

i-kk 



(en) 



% (i + Wl ^)(fc2 + e ^ L ) 
--2 



t>-2 
?POL 



1 + 



POL 



where the Fourier integral is given by @j and Cpol ^ s ^ ne 
effective screening length in the polymer solution. 

Averaging over the distribution of the segments and 
integrating over the solid angle, recalling that W POL (fc) 
does not depend on the direction of k , we get : 

1 



3tH J 



k 1 



770/c 2 ~W[ OL (k) 



x exp 



k 2 l\s 



(> 



(C12) 



Then, applying the Fourier transform defined at the be- 
ginning of the appendix, one finds : 



1 



dk 



k 2 



V(fc)i 

(C13) 



where we use d fd flt (1 — kk) = ^1 and l qq > has been 
evaluated in (C1C). 

In order to make the calculation of Q analytically 
tractable, without coarsening the physical picture too 
much, we will make two approximations : 

1. the Kirkwood-Riseman approximation, which 
amounts to discard the off-diagonal terms in T qq , (k) . 

2. the long chain limit, n^> 1. 

Under these assumptions, both Q and its inverse are di- 
agonal and their Fourier elements are simply related as : 



{g- l ) qq > = {Q- 1 ) q 5 qql = {g q y 1 5 qq , 



(C14) 



For convenience we will denote these matrix elements by 

We will proceed further to evaluate W POL (k) . Notic- 
ing that I q = I- q , one finds : 



W 



POL, 



(k)^-c POL \\ {k)g^ + 2Y J l q (k)G q 1 ^ 1 (C15a) 
W POi (fc) = c POL Io(fe) Gv 1 1 (C15b) 
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in which we obtained (g t from 



ft) = E ( K ^) = j-2 [ L f L d s d s 'g-\ s , s ') = g^ 



L 2 



o Jo 



<ft X > = 0o 



(C16) 



Working in the hydrodynamic limit k — > 0, we expand 
l q (k <C 1) — -^zk 2 R 2 and converting the sum over q 

to an integral over the wavevector ji = 2 -p, we obtain 
the foll owing expre ssion for W POL (k) by adding the two 
terms ( |cT5a| , pl5h| ) : 



W 



POL 



,OL R 2 g k 2 



dfi — g^t (ci7) 



1 



ax 



3tt 2 770(1 + W POL ) R g (3 J x 2 + 1 
1 



x <^ 1 



{fix) 



[1 - exp(-/? 2 ^)] ^ 1 (C22) 



with /3 a dimensionless variable : (3 = RgS,' pol- The 
integral gives the exact result : 

1 1 1 *(fi) + 



3tt m (l + Wf OL ) R g (i 



P(/3) = l + ^ 



1 - exp(/3 2 ) + -= exp(/3 2 )r(|, 0, /3 2 ) 



(C23) 



r(|,0,a;) = / dtt* exp(-i) 



To calculate the translational friction coefRcient (and 
implicitlythe diffusion coefficient), we have to consider 
one polymer chain moving with some uniform velocity 
Uo inside an effective medium replacing all others Np — f 
chains, medium char acterized by the force propagator 
£?(k) given by (Cll). The equation of motion of the 
center of mass of the chain is : 



MR 



(C18) 



where the last expression is the incomplete Gamma func- 
tion. 

Thus the friction coefficient is given by : 

C* i = Go 1 = tern [i + wf UL ($ 



(C24) 



From ( |C13| ), in the long-chain limit, Q ^ becomes 



M is the mass of the chain and er^ is the force exerted by 
the i-th bead of the polymer upon the fluid. Neglecting 
the inertial term (stationary problem) in the above equa- 
tion of motion, the total external average force upon the 
polymer chain is found to be : 



(fo) = (E-0=E( /c ^>-uo 



(C19) 



where the expression of cr.j is a particular case of (|A4| ) 
no spheres present and the other chains are absorbed in 
the effective fluid by replacing K" 1 with JCy . The trans- 
lational friction coefficient follows readily : 



(C20) 



where first we had a configurational average as defined 
in ([32]) that reduced to an average over the distribution 
of the segments S^, about the center of mass of the 
chain. Introducing the arc-length variables s and s' and 
applying the Fourier transform we obtain for £ t : 



Go 1 1 - (g t ) 1 



(C21) 



Next we insert (CIO) in ( C13j ) and setting q = q' = we 
calculate Qq : 



1 



I 



To. \ d P- 



P 



-2 

POL 



3ir 2 Vo (l + W POL ) L 2 J a V+£' 

(I (C25) 



1 p 2 R 2 9 



where \x = -j^. Working out the integral in terms of 
dimensionless variables we get : 



V2 



V2- 



3/2 



6tt 770(1 + W 1 POL ) R g [3 



(C26) 



Also, Q 1 is immediately computed as : Q = (Q M )~ 



Finally, changing the integration variable in (C17) 
from (j, to x = (fiL/f3 2 ) and defining the volume fraction 

of the polymers to be <& POL = 47r ^ 3 c POL , we calculate 
W POi (k) : 



W POL (k) = -- T70 (1 + W POL ) <S> POL k 2 \ 

7T p 



d.v 



x 4 + l 



V2x 3 (x 3 -X + V2) 



t (C27) 
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where the integration can be done exactly, leading even- 
tually to the desired equation for the self-energy W POL 
of a suspension of Np non- interacting, mobile polymer 
chains : 



W 



POL 



(k) = --no (1 + W[ OL ) $ POL k 2 Q(J3) 1 (C28) 



with (5 = 



As a last note, we present here the series expansions 
and the asymptotic behavior of the functions P(x) and 
Q(x) introduced in the previous calculations : 



x < 1 



P(x) 



Q(x) ~ 2 



-2- -4 



l x 2 + 8 a 

2 X T 15v^ 



1 ^ _ 1 ^3 + 

3-n-\/-7r 47T- 1 



X > 1 



P( X ) ~ 1- 2 +J? _ 1, 



Q(x) ~ ±x + 



(C29) 
(C30) 



+ ... 
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FIG. 1. The physical interpretation of the multiple scat- 
tering expression of the microscopic velocity field v(r). O = 
observation point; G = the Oseen tensor of the pure fluid; 
T = single-object velocity propagator. 



FIG. 11. Differences in the diffusion coefficient of the poly- 
mer chain function of $ 5P for four pairs of t values (1P/S 
system). From top to bottom,following the maximum of each 
curve : O D(t = 0.01) - D(t = 10); O D(t = 1) - D(t = 10); 
A D(t = 0.01) - D(t = 1); □ D(t = 0.01) - D(t = 0.1). 



FIG. 2. Diffusion coefficient D function of $^ VL for the 
1S/P limit case and for the S/P system with $ SP param- 
eter; all t values. From top to bottom : + 1 sphere limit; 
O $ SP = 0.01; A $ SP = 0.11; O $ SP = 0.21; □ $ SP = 0.31. 



FIG. 3. Diffusion coefficient D function of <J> for a probe 

sphere moving among other fixed spheres in a suspension of 

polymers (S/P system), with $ POL parameter; all t values. 

From top to bottom : O $ POL = 0.01; A $ POi = 0.09; 
q $pol = 16 . n qPOl = 24 



FIG. 12. Diffusion coefficient D function of $ POi for one 
mobile polymer chain inside a dispersion of other polymer 
chains and fixed spheres (P/S system), with 3> SP parameter; 
t = 0.1. From top to bottom , Figs. [12-14] : O $ 5P = 0.01; 
A $ SP = 0.14; O $ SP = 0-26; □ $ 5P = 0.39. 



FIG. 13. Diffusion coefficient D function of $> POL . P/S 
system, with $ SP parameter; t = 1. 



FIG. 14. Diffusion coefficient D function of $ POi . p/S 



FIG. 4. Total viscosity change 2— 
for the 1S/P limit case and for the S/P system with $ 5P pa- 
rameter; all t values. From bottom to top : + 1 sphere limit; 

0.26; □ $ SP = 0.39. 



1 function of $ POL system, with $ parameter; t = 10. 



O $ SP = 0.01; A $ SP = 0.14; Q $ £ 



FIG. 5. Total viscosity change 1— 

for the S/P system, with $ POi 
O $ POi : 



parameter; all t values. 



From bottom to top 
O 4> POL =0.16; □ <S> POL 



0.01; A 



0.09; 



0.24. 



FIG. 6. Relative viscosity change 



1 function of 



<& POL for the S/P system, with $ SP parameter; all t val- 
ues. From bottom to top : O $ SP = 0.09; A $ 5P = 0.24; 
O $ SP = 0.39; □ $ 5P = 0.46. 



FIG. 7. Diffusion coefficient D function of $ 5P for the 

1P/S limit and for the P/S system with $ POL parameter 

t — 0.01. From top curve to bottom curve, Figs. [7-10] 

+ 1 chain limit; O $ POi = 0.01; A $ POL = 0.09 
q $pol = 16 . n ^POL = 24 



FIG. 15. Total viscosity change ^- - 1 function of $ POi . 
P/S system, with $ SP parameter; t = 0.01. From bottom to 
top curve : O $ 5P = 0.01; A $ SP = 0.14; Q ® SP = 0.26; 



1 function of $ SP □ 



0.39. 



FIG. 16. Total viscosity change 



n 



P/S system, with <E> parameter; t 

top curve : O $ 5P = 0.01; A $ SP = 0.14; Q = 0.26 



- 1 function of $ POi . 
10. From bottom to 

SP 



□ $- 



0.39. 



FIG. 17. Total viscosity change 2!!! - 1 function of $ SP . 
P/S system, with $ POi parameter; t = 0.01. From bottom 

: 0.01; A $ POi = 0.09; 



to top curve, Figs. [17-18] 
O $ POL =0.16; □ $ POi 



O $ POL 
= 0.24. 



eff q p 

FIG. 18. Total viscosity change 2— - 1 function of $ 6 . 
P/S system, with $ POi parameter; t = 10. 



FIG. 8. Diffusion coefficient D function of $ SP . 1P/S 
and P/S systems, with $ POL parameter; t = 0.1. 



FIG. 9. Diffusion coefficient D function of $ SP . 1P/S 
and P/S systems, with $ POL parameter; t = 1. 



FIG. 19. The relative viscosity change 3 1 as a func- 
tion of $ POi for the P/S system, with $ SP parameter; t = 1 
when symbols are joined with lines and t = 10 when they 
are isolated. The corresponding $ SP values are : O and 

O-O $ SP = 0.01; A $ SP = 0.14; O * 5P = 0.26; □ and 
□_□ $sp = q 39 



FIG. 10. Diffusion coefficient D function of $ 5P . 1P/S 
and P/S systems, with $ POL parameter; t = 10. 



FIG. 20. The relative viscosity change — 1 as a function 
of <E> SP . P/S system, with $ POL parameter; t = 0.01. From 

bottom to top curve, A $ POL = 0.09; O * POL = 0.16; 
D §pol = 24 



24 



eff 

FIG. 21. The relative viscosity change 3—^ — 1 as a func- 
tion of $ SP . P/S system, with $ POi parameter; t = 1 
when symbols are joined with lines and £ = 10 when they 
are isolated. The corresponding <E> P values are : A and 
A-A $ POL = 0.09; O and OO $ POL = 0.16; □ and 
□_□ $ SP = 0.24. 
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